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PREFACE 

This  report  was  prepared  by  Dr.  George  D.  Ashton,  Chief,  Geophysical  Sciences  Branch,  Re¬ 
search  Division.  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory.  Sandra  J.  Smith 
processed  much  of  the  data.  Anthony  Gow  and  John  Govoni  provided  ready  access  to  the  original, 
detailed  thickness  data.  The  report  was  technically  reviewed  by  Michael  Bilello  and  Darryl  Cal¬ 
kins.  The  work  was  funded  by  project  CWIS  3 1 750,  Prediction  of  Ice  Formation. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional  purposes.  Citation 
of  brand  names  does  not  constitute  an  official  endorsement  or  approval  of  such  commercial 
products. 


PREDICTING  LAKE  ICE  DECAY 


George  D.  Ashton 


BACKGROUND 

In  dealing  with  the  effects  of  ice  on  water  bodies 
or  in  assessing  ice  problems,  the  most  important 
characteristic  of  the  ice  cover  is  its  thickness.  Over 
the  years  much  attention  has  been  devoted  to  pre¬ 
dicting  the  thickening  of  the  ice  cover  under  the  in¬ 
fluence  of  subfreezing  air  temperatures,  but  relatively 
little  attention  has  been  devoted  to  the  decay  of 
thickness  during  the  spring.  In  part  the  greater  at¬ 
tention  to  the  thickening  stage  is  due  to  the  avail¬ 
ability  of  a  simple  heat-conduction  model  by  which 
to  predict  the  growth  rates  (the  “Stefan”  problem); 
a  simple  analysis  yields  the  well-known  “square-root- 
of-freezing-degree-days”  method.  During  the  decay 
period  both  the  top  and  the  bottom  of  a  floating 
ice  cover  are  at  0°C;  hence,  the  heat-conduction 
framework  of  analysis  is  not  applicable,  and  the 
problem  of  predicting  thinning  rates  is  conceptually 
unclear.  This  report  presents  a  simple  algorithm  for 
predicting  decay  of  lake  ice  and  tests  it  on  data  from 
a  nine-year  study  of  the  lake  ice  cover  on  Post  Pond 
in  Lyme,  New  Hampshire. 

There  is  little  published  work  on  the  decay  per¬ 
iod  of  lake  ice  covers.  Bileilo  (1980)  presented  ex¬ 
tensive  data  from  a  variety  of  locations  and  examined 
several  empirical  approaches  to  predict  decay.  He 
found  that  an  accumulated  thawing  degree-day 
(ATDD)  index  provided  the  best  overall  correlations, 
particularly  when  the  ice  cover  was  not  subjected  to 
water  currents  or  other  action  that  would  mechan¬ 
ically  break  up  the  ice  cover.  The  relationship  used 
by  Bileilo  (1980)  is 


where 

/  =  thickness  of  the  decaying  ice  sheet 
(cm) 

/  =  maximum  ice  thickness  at  the  start 
of  ice  decay  (cm) 

S  =  “slope”  of  the  ice  thickness  vs  ATDD 
index  plot  (cm  0C~* ) 

ATDD  =  accumulated  thawing  degree-days 
(°C  with  0°C  base). 

This  report  uses  a  similar  approach  but  puts  it  in  a 
more  rigorous  form  by  performing  a  simple  heat 
balance  to  obtain  the  same  result. 


ANALYSIS 

This  analysis  consists  of  equating  the  energy  flux 
associated  with  a  melting  rate  of  t/rj/dr  to  a  sensible 
heat  flux  relationship  based  on  the  difference  be¬ 
tween  0°C  and  the  above-freezing  air  temperature; 
that  is, 

(2) 

where 

pj  =  density  of  the  ice  (916  kg  m~3) 

A  =  heat  of  fusion  (3.34  x  10s  J  kg'1) 
dntdt  =  thinning  rate  (m  s'1 )  of  an  ice  cover 
of  thickness  r? 
t  =  time 

Hti  =  transfer  coefficient  (W  m~2  °C~l ) 

Tm  =  melting  temperature  (0°C) 

T,  =  air  temperature. 
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Integration  of  eq  2  yields 


DISCUSSION  OF  RESULTS 
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where  C  is  the  constant  of  integration.  If  tjm  is  the 
ice  thickness  at  the  beginning  of  melt,  then  C  =  r?m, 
the  term  (7^  -  Ta)  t  is  the  ATDD  (with  appropriate 
units),  and  eq  3  is  equivalent  to  eq  1 . 

For  lake  ice  in  Canada  and  Alaska,  Bilello  found 
values  of  5 in  eq  1  (or  Hai /PjA  in  eq  3)  to  vary  be¬ 
tween  0.43  x  10"7  m  s"1  “C"1  (0.37  cm  day-1  “C"1 ) 
and  1 .49  x  10"7  m  s"1  °C"1  (1 .29  cm  day"1  °C_1 ), 
with  an  average  value  for  IS  lakes  of  0.79  x  10'7 
ms"1  “C"1  (0.68  cm  day"1  “C"1).  Not  surprisingly 
the  equivalent  values  of  //ai  are  very  similar  to  typ¬ 
ical  heat-transfer  coefficients  for  sensible  heat  trans¬ 
fer  from  air  to  water  surfaces,  which  are  on  the  order 
of  15-25  Wm"2  °C-'. 


Examination  of  ail  nine  cases  leads  to  the  con¬ 
clusion  that  a  reasonable  value  for  Hai  for  this  lake 
is  between  1 5  and  20  W  m"2  “C"1 .  The  use  of  equa¬ 
tion  5  for  cold  spells  (<  0°C)  seems  to  adequately 
treat  delays  of  thawing  or  slight  increases  in  thick¬ 
ness.  Overall,  the  agreement  between  calculated  and 
observed  thicknesses  is  adequate,  particularly  in  light 
of  the  extreme  simplicity  of  the  method.  To  include 
the  processes  and  detailed  calculations  necessary  to 
improve  the  method  would  require  that  the  calcula¬ 
tions  be  an  order  of  magnitude  more  complex  and 
that  the  input  data,  such  as  wind  speed,  radiation 
components  and  cloudiness  be  of  much  greater  de¬ 
tail.  For  many  applications  the  simple  method  pre¬ 
sented  here  should  prove  adequate  and  (for  short¬ 
term  predictions)  much  better  than  historical  aver¬ 
ages. 


APPLICATION  TO  POST  POND  DATA 

Post  Pond  is  a  small  eutrophic  lake  of  glacial 
origin  located  in  Lyme,  New  Hampshire,  at  a  lati¬ 
tude  of  44°.  Its  area  is  0.46  km2  and  its  maximum 
depth  is  12  m.  The  ice  thickness  measurements 
used  here  were  obtained  during  the  winters  from 
1973-74  to  1981-82  and  are  reported  more  fully 
by  Gow  and  Govoni  (1983).  Ice  thickness  was 
measured  about  once  a  week. 

The  present  method  begins  with  the  maximum 
ice  thickness  of  each  season  and  calculates  the  ex¬ 
pected  daily  thickness  change  using  eq  2  written  in 
the  form 

Atj  =  <rm  -  TJ  *  <4> 

where  Ta  was  taken  as  the  average  of  the  daily  max¬ 
imum  and  minimum  air  temperatures  measured  at 
Hanover,  New  Hampshire,  some  20  km  south  of 
Post  Pond.  Figure  1  provides  a  comparison  of  pre¬ 
dictions  with  observations  for  Hti  values  of  10,  15 
and  20  W  m"2  “C"1 .  When  the  ice  was  thickening 
(i.e.  when  Tt  <  0°C),  the  algorithm  used  was  (Ash¬ 
ton  1980) 


LIMITATIONS 

The  extreme  simplicity  of  this  method  of  pre¬ 
dicting  lake  ice  decay  necessarily  has  limitations. 

One  limitation  is  that  the  appropriate  value  for  the 
transfer  coefficient  has  been  determined  for  one 
site  only  and  may  not  be  appropriate  for  sites,  times, 
or  locations  where  other  components  of  the  energy 
budget  are  of  different  relative  values.  The  method 
clearly  does  not  apply  near  the  shorelines  or  where 
stream  inlet  flows  cause  enhanced  melting. 

Implicit  in  the  method  is  the  assumption  that  the 
ice  surface  is  reasonably  well  drained  of  meltwater. 
This  is  usually  true,  either  as  a  consequence  of  the 
porosity  that  develops  inside  the  ice  cover  during 
the  deterioration  process  or  as  a  result  of  natural 
cracks  or  holes  in  the  ice  cover.  Wake  and  Rumer 
( 1 979)  discussed  some  of  the  difficulties  that  occur 
when  the  ice  is  not  well  drained. 

This  method  does  not  consider  the  snow  on  top 
of  the  ice  cover.  During  the  initial  melt  period  and 
after  subsequent  snowfalls  this  may  lead  to  some  in¬ 
accuracy,  although  generally  during  spring  the  snow 
(but  not  the  snow  ice)  rapidly  melts  to  the  top  of 
the  ice. 

The  method  also  should  not  be  applied  for  per¬ 
iods  less  than  a  day  or  so  because  of  the  highly  var¬ 
iable  effects  of  radiation  exchange  that  occur  on  a 
diurnal  basis.  Finally  predictions  based  on  this 
method  are  only  as  good  as  the  forecasts  of  air 
temperature. 
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Figure  1.  Predicted  ice  thicknesses  during  the  decay  period  starting  from  maximum  observed  ice  cover 
thickness  and  using  different  Ho(  values.  The  plotted  points  are  observed  thicknesses. 
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Figure  1  (cant'd).  Predicted  ice  thicknesses  during  the  decay  period  starting  from  maximum  observed  ice 
cover  thickness  and  using  different  H4(  values.  The  plotted  points  are  observed  thicknesses. 
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